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Abstract 

The considerable part of the parameter space in the MSSM corresponding to the 
infrared quasi fixed point scenario is excluded by LEPII bounds on the lightest 
Higgs boson mass. In the NMSSM the mass of the lightest Higgs boson reaches its 
maximum value in the strong Yukawa coupling limit when Yukawa couplings are 
essentially larger than gauge ones at the Grand Unification scale. In this case the 
renormalization group flow of Yukawa couplings and soft SUSY breaking terms is 
investigated. The quasi-fixed and invariant lines and surfaces are briefly discussed. 
The coordinates of the quasi-fixed points, where all solutions are concentrated, are 
given. 
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1 Introduction 



The search for the Higgs boson remains one of the top priorities for existing accelerators 
as well as for those still at the design stage. This is because this boson plays a key role 
in the Standard Model which describes all currently available experimental data with a 
high degree of accuracy. As a result of the spontaneous symmetry breaking SU (2) ®U{1) 
the Higgs scalar acquires a nonzero vacuum expectation value without destroying the 
Lorentz invariance, and generates the masses of all fermions and vector bosons. An 
analysis of the experimental data using the Standard Model has shown that there is a 
95% probability that its mass will not exceed 210 GeV jlj. At the same time, assuming 
that there are no new fields and interactions and also no Landau pole in the solution 
of the renormalisation group equations for the self-action constant of Higgs fields up to 
the scale Mpi ^ 2.4 ■ 10^*^ GeV, we can show that mh < 180 GeV 0,0. In this case, 
physical vacuum is only stable provided that the mass of the Higgs boson is greater than 
135 GeV [21- El- However, it should be noted that this simplified model does not lead to 
unification of the gauge constants JT^ and a solution of the hierarchy problem [HI. As a 
result, the construction of a realistic theory which combines all the fields and interactions 
is extremely difficult in this case. 

Unification of the gauge constants occurs naturally on the scale Mx ~ 3 ■ 10^^ GeV within 
the super symmetric generalisation of the Standard Model, i.e., the Minimal Supersym- 
metric Standard Model (MSSM) 0. In order that all the fundamental fermions acquire 
mass in the MSSM, not one but two Higgs doublets Hi and H2 must be introduced 
in the theory, each acquiring the nonzero vacuum expectation value Vi and V2 where 
v"^ = vl + V2 = (246 GeV)^. The spectrum of the Higgs sector of the MSSM contains 
four massive states: two CP-even, one CP-odd, and one charged. An important distin- 
guishing feature of the supersymmetric model is the existing of a light Higgs boson in 
the CP-even sector. The upper bound on its mass is determined to a considerable extent 
by the value tan/3 = v^jvx. In the tree- level approximation the mass of the lightest 
Higgs boson in the MSSM does not exceed the mass of the Z-boson [Mz ~ 91.2 GeV): 
"^h < Mz \ cos2/3| jH]. Allowance for the contribution of loop corrections to the effective 
interaction potential of the Higgs fields from a t-quark and its superpartners significantly 
raises the upper bound on its mass: 



Here A are the loop corrections JO], JT]. The values of these corrections are proportional 
to mj, where nit is the running mass of t-quark which depends logarithmically on the 
supersymmetry breaking scale Ms and is almost independent of the choice of tan/3. 
In bounds on the mass of the Higgs boson were compared in the Minimal 

Standard and Supersymmetric models. The upper bound on the mass of the light CP- 
even Higgs boson in the MSSM increases with increasing tan/3 and for tan/3 ^ 1 in 
realistic supersymmetric models with Ms < 1000 GeV reaches 125 — 128 GeV. 

However, a considerable fraction of the solutions of the system of MSSM renormalisation 
group equations is focused near the infrared quasi-fixed point at tan /? ~ 1. In the region 
of parameter space of interest to us (tan/? -C 50) the Yukawa constants of a 6-quark 
(/lb) and a r-lepton [h^) are negligible so that an exact analytic solution can be obtained 
for the one-loop renormalisation group equations For the Yukawa constants of a 
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t-quark ht{t) and the gauge constants gi{t) its solution has the following form: 
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F{t) = j E{t')dt' 



where the index i has values between 1 and 3, 



hi = 33/5, 62 = 1, h = -3 
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The variable t is determined by a standard method t = ln(M|-/g^). The boundary 
conditions for the renormalisation group equations are usually set at the grand unifi- 
cation scale Mx {t = 0) where the values of all three Yukawa constants are the same: 
ai{0) = a2(0) = 153(0) = ctQ. On the electroweak scale where hf{0) ^ 1 the second term 
in the denominator of the expression describing the evolution of Yt{t) is much smaller 
than unity and all the solutions are concentrated in a narrow interval near the quasi-fixed 
point yQFp(t) = E{t)/6F{t) [TH]. In other words in the low-energy range the dependence 
of Yt{t) on the initial conditions on the scale Mx disappears. In addition to the Yukawa 
constant of the t-quark, the corresponding trilinear interaction constant of the scalar 
fields At and the combination of the scalar masses dJlf = rriQ + mfj + rn^ also cease to 
depend on At{0) and dJlf{0) as Yt{0) increases. Then on the electroweak scale near the 
infrared quasi-fixed point At{t) and OJt^(t) are only expressed in terms of the gaugino 
mass on the Grand Unification scale. Formally this type of solution can be obtained if 
Yt{0) is made to go to infinity. Deviations from this solution are determined by ratio 
l/6F{t)Yt{0) which is of the order of l/10hf{0) on the electroweak scale. 

The properties of the solutions of the system of MSSM renormalisation group equations 
and also the particle spectrum near the infrared quasi-fixed point for tan /5 ~ 1 have 
been studied by many authors Recent investigations have shown that 

for solutions Yt{t) corresponding to the quasi- fixed point regime the value of tan/5 is 
between 1.3 and 1.8. These comparatively low values of tan/3 yield significantly more 
stringent bounds on the mass of the lightest Higgs boson. The weak dependence of the 
soft supersymmetry breaking parameters At{t) and OJlj(t) on the boundary conditions 
near the quasi-fixed point means that the upper bound on its mass can be calculated 
fairly accurately. A theoretical analysis made in [THj , [lEl showed that ruh does not exceed 
94 ± 5 GeV. This bound is 25 — 30 GeV below the absolute upper bound in the Minimal 
Supersymmetric Model. Since the lower bound on the mass of the Higgs boson from 
LEPII data is 113 GeV [Tj, which for the spectrum of heavy supersymmetric particles is 
the same as the corresponding bound on the mass of the Higgs boson in the Standard 
Model, a considerable fraction of the solutions which come out to a quasi-fixed point 
in the MSSM, are almost eliminated by existing experimental data. This provides the 
stimulus for theoretical analyses of the Higgs sector in more complex supersymmetric 
models. 
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The simplest expansion of the MSSM which can conserve the unification of the gauge 
constants and raise the upper bound on the mass of the hghtest Higgs boson is the 
Next-to-Minimal Supersymmetric Standard Model (NMSSM) [IH|-|2ni- By definition the 
superpotential of the NMSSM is invariant with respect to the discrete transformations 
y'a — ^'^^^'^Va of the Z3 group [T^l which means that we can avoid the problem of the 
/^-term in supergravity models. The thing is that the fundamental parameter /i should 
be of the order of Mpi since this scale is the only dimensional parameter characterising 
the hidden (gravity) sector of the theory. In this case, however, the Higgs bosons H\ 
and Hi acquire an enormous mass rn\^^ ^ ^ ^^"^ breaking of SU{2) Cg) U{1) 
symmetry occurs. In the NMSSM the term ^{HiHi) in the superpotential is not invariant 
with respect to discrete transformations of the Z3 group and for this reason should be 
eliminated from the analysis (/x = 0). As a result of the multiplicative nature of the 
renormalisation of this parameter, the term fi{q) remains zero on any scale q < MxH-Mpi. 
However, the absence of mixing of the Higgs doublets on electroweak scale has the result 
that Hi acquires no vacuum expectation value as a result of the spontaneous symmetry 
breaking and d-type quarks and charged leptons remain massless. In order to ensure that 
all quarks and charged leptons acquire nonzero masses, an additional singlet superfield Y 
with respect to gauge SU{2) ® U{1) transformations is introduced in the NMSSM. The 
superpotential of the Higgs sector of the Nonminimal Supersymmetric Model [THj-pU] 
has the following form: 

Wj, = \Y{HiHi) + (3) 

As a result of the spontaneous breaking of SU{2) ® U{1) symmetry, the field Y acquires 
a vacuum expectation value ((F) = y/V^.) and the effective /i-term (/x = Xy/\^) is 
generated. 

In addition to the Yukawa constants A and x, and also the Standard Model constants, 
the Nonminimal Supersymmetric Model contains a large number of unknown parame- 
ters. These are the so-called soft supersymmetry breaking parameters which are required 
to obtain an acceptable spectrum of superpartners of observable particles form the phe- 
nomenological point of view. The hypothesis on the universal nature of these constants 
on the Grand Unification scale allows us to reduce their number in the NMSSM to three: 
the mass of all the scalar particles mo, the gaugino mass M1/2, and the trilinear inter- 
action constant of the scalar fields A. In order to avoid strong CP-violation and also 
spontaneous breaking of gauge symmetry at high energies (Mpi ^ E ^ nit) as a result 
of which the scalar superpartners of leptons and quarks would require nonzero vacuum 
expectation values, the complex phases of the soft supersymmetry breaking parameters 
are assumed to be zero and only positive values of mg are considered. Naturally universal 
supersymmetry breaking parameters appear in the minimal supergravity model ^S] and 
also in various string models j22],j2l]. In the low-energy region the hypothesis of universal 
fundamental parameters allows to avoid the appearance of neutral currents with fiavour 
changes and can simplify the analysis of the particle spectrum as far as possible. The 
fundamental parameters thus determined on the Grand Unification scale should be con- 
sidered as boundary conditions for the system of renormalisation group equations which 
describes the evolution of these constants as far as the electroweak scale or the super- 
symmetry breaking scale. The complete system of the renormalisation group equations 
of the Nonminimal Supersymmetric Model can be found in [22] , |2E] • These experimental 
data impose various constraints on the NMSSM parameter space which were analysed in 

E71,|2H1. 
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The introduction of the neutral field Y in the NMSSM potential leads to the appearance of 
a corresponding F-term in the interaction potential of the Higgs fields. As a consequence, 
the upper bound on the mass of the lightest Higgs boson is increased: 



mh<y —v"^ sin^ 2/5 + M| cos^ 2(3 + A^V + ■ (4) 

The relationship Q was obtained in the tree- level approximation (An = 0) in pUj. 
However, loop corrections to the effective interaction potential of the Higgs fields from 
the t-quark and its superpartners play a very significant role. In terms of absolute 
value their contribution to the upper bound on the mass of the Higgs boson remains 
approximately the same as in the Minimal Supersymmetric Model. When calculating 
the corrections A^^-* and A^^-* we need to replace the parameter /i by \y/\/2. Studies of 
the Higgs sector in the Nonminimal Supersymmetric model and the one-loop corrections 
to it were reported in |2II,[2S],|2S1-[SI1- In [Hj the upper bound on the mass of the 
lightest Higgs boson in the NMSSM was compared with the corresponding bounds on rrih 
in the Minimal Standard and Supersymmetric Models. The possibility of a spontaneous 
CP-violation in the Higgs sector of the NMSSM was studied in 

It follows from condition that the upper bound on irih increases as A increases. More- 
over, it only differs substantially from the corresponding bound in the MSSM in the range 
of small tan (3. For high values (tan /5 S> 1) the value of sin 2/? tends to zero and the upper 
bounds on the mass of the lightest Higgs boson in the MSSM and NMSSM are almost 
the same. The case of small tan j3 is only achieved for fairly high values of the Yukawa 
constant of a t-quark ht on the electroweak scale {ht{to) > 1 where to = ln(M|^/m|)), and 
tan (3 decreases with increasing ht{tQ). However, an analysis of the renormalisation group 
equations in the NMSSM shows that an increase of the Yukawa constants on the elec- 
troweak scale is accompanied by an increase of /?.t(0) and A(0) on the Grand Unification 
scale. It thus becomes obvious that the upper bound on the mass of the lightest Higgs 
boson in the Nonminimal Supersymmetric model reaches its maximum on the strong 
Yukawa coupling limit, i.e., when /it(0) ^ fl'j(O) and A(0) ^ fi'j(O). 



2 Renormalization of the Yukawa couplings 

From the point of view of a renormalisation group analysis, investigation of the NMSSM 
presents a much more complicated problem than investigation of the minimal SUSY 
model. The full set of renormalization group equations within the NMSSM can be found 
in |2SI,[2n]- Even in the one-loop approximation, this set of equations is nonlinear and 
its analytic solution does not exist. All equations forming this set can be partitioned into 
two groups. The first one contains equations that describe the evolution of gauge and 
Yukawa coupling constants, while the second one includes equations for the parameters 
of a soft breakdown of SUSY, which are necessary for obtaining a phenomenologically 
acceptable spectrum of superpartners of observable particles. Since boundary conditions 
for three Yukawa coupling constants are unknown, it is very difficult to perform a nu- 
merical analysis of the equations belonging to the first group and of the full set of the 
equations given above. In the regime of strong Yukawa coupling, however, solutions to 
the renormalisation group equations are concentrated in a narrow region of the parameter 
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space near the electroweak scale, and this considerably simplifies the analysis of the set 
of equations being considered. 

In analysing the nonlinear differential equations entering into the first group, it is conve- 
nient to use the quantities pt, p\, p^, pi, and p2, defined as follows: 

''^'^-^y ''^'^-^y '-^'^-^y ''^'^'^y ''^'^'^y 

where a,{t) = ^7K^)/(4v^)^ Y,{t) = /i?(t)/(47r)2, Y,{t) = X\t)/{A7ry, and 
Y^it) = x2(t)/(47r)2. 

Let us first consider the simplest case of x = 0. The growth of the Yukawa coupling 
constant A(to) at a fixed value of /it (to) results in that the Landau pole in solutions to 
the renormalization group equations approaches the Grand Unification scale from above. 
At a specific value A(to) = Amax, perturbation theory at g ~ Mx cease to be applicable. 
With increasing (decreasing) Yukawa coupling constant for the 6-quark, Amax decreases 
(increases). In the {pt, P\) plane, the dependence Xl^g^{h^) is represented by a curve 
bounding the region of admissible values of the parameters pt(to) and pxitQ). At px = 0, 
this curve intersects the abscissa at the point pt = p^^^(to)- This is the way in which 
there arises, in the {pt,px) plane, the quasi-fixed (or Hill) line near which solutions to 
the renormalization group equations are grouped (see Fig. 1). With increasing A^(0) and 
hf{0), the region where the solutions in questions are concentrated sharply shrinks, and 
for rather large initial values of the Yukawa coupling constants they are grouped in a 
narrow stripe near the straight line 

Pt(to) + 0.506pA(to) = 0.91, (5) 

which can be obtained by fitting the results of numerical calculations (these results are 
presented in Fig. 2). Moreover, the combination /i^(to) + 0.506A^(to) of the Yukawa 
coupling constants depends much more weakly on A^(0) and hf{0) than A^(to) and hlito) 
individually IHE]- In other words, a decrease in A^(to) compensates for an increase in 
hflto), and vice versa. The results in Fig. 3, which illustrate the evolution of the above 
combinations of the Yukawa coupling constants, also confirm that this combination is 
virtually independent of the initial conditions. 

In analysing the results of numerical calculations, our attention is engaged by a pro- 
nounced nonuniformity in the distribution of solutions to the renormalization group 
equations along the infrared quasi-fixed line. The main reason for this is that, in the 
regime of strong Yukawa coupling, the solutions in question are attracted not only to 
the quasi-fixed but also to the infrared fixed (or invariant) line. The latter connects 
two fixed points. Of these, one is an infrared fixed point of the set of renormalization 
group equations within the NMSSM {pt = 7/18, px = 0, pi = 0, p2 = 0) [SB], while the 
other fixed point {px/ Pt = 1) corresponds to values of the Yukawa coupling constants in 
the region Yt,Yx 3> ai, in which case the gauge coupling constants on the right-hand 
sides of the renormalization group equations can be disregarded [34J. For the asymptotic 
behaviour of the infrared fixed line at pt, pa ^ 1 we have 

8 2 

^^ = ^'-15-75^^' 
while in the vicinity of the point pt = 7/18, pa = we have 

Px ~ (pt - 7/18)25/14. 
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The infrared fixed line is invariant under renormalization group transformations - that is, 
it is independent of the scale at which the boundary values Yt{0) and 1a(0) are specified 
and of the boundary values themselves. If the boundary conditions are such that Yt{0) and 
Yx{0) belong to the fixed line, the evolution of the Yukawa coupling constants proceeds 
further along this line toward the infrared fixed point of the set of renormalization group 
equations within the NMSSM. With increasing t, all other solutions to the renormalization 
group equations are attracted to the infrared fixed line and, for t/(47r) ^ 1, approach the 
stable infrared fixed point. From the data in Figs. 1 and 2, it follows that, with increasing 
Yt{0) and Yx{0), all solutions to the renormalization group equations are concentrated in 
the vicinity of the point of intersection of the infrared fixed and the quasi-fixed line: 

pf^'ito) = 0.803, pf P(to) = 0.224. 

Hence, this point can be considered as the quasi-fixed point of the set of renormalization 
group equations within the NMSSM at x = 0. 

In a more complicated case where all three Yukawa coupling constants in the NMSSM 
are nonzero, analysis of the set of renormalization group equations presents a much more 
difficult problem. In particular, invariant (infrared fixed) and Hill surfaces come to the 
fore instead of the infrared fixed and quasi-fixed lines. For each fixed set of values of 
the coupling constants Yt{to) and Y^(to), an upper limit on Yx(^o) can be obtained from 
the requirement that perturbation theory be applicable up to the Grand Unification scale 
Mx- A change in the values of the Yukawa coupling constants ht and x at the electroweak 
scale leads to a growth or a reduction of the upper limit on 1a (to)- The resulting surface 
in the {pt,p^,p\) space is shown in Fig. 4. In the regime of strong Yukawa coupling, 
solutions to the renormalization group equations are concentrated near this surface. In 
just the same way as in the case of = 0, a specific linear combination of Yt, Y\, and 
Y^ is virtually independent of the initial conditions for Yi{0) oo: 

Pt{to) + 0.72px{to) + 0.33p^(to) = 0.98. (6) 

The evolution of this combination of Yukawa couplings at various initial values of the 
Yukawa coupling constants is illustrated in Fig. 5. 

On the Hill surface, the region that is depicted in Fig. 4 and near which the solutions in 
question are grouped shrinks in one direction with increasing initial values of the Yukawa 
coupling constants, with the result that, at Yt{0), Y^{0), and 1a(0) ~ 1, all solutions 
are grouped around the line that appears as the result of intersection of the quasi-fixed 
surface and the infrared fixed surface, which includes the invariant lines lying in the 
p^ = and Pa = planes and connecting the stable infrared point with, respectively, 
the fixed point px/ pt = 1 and the fixed point p^/pt = 1 in the regime of strong Yukawa 
coupling. In the limit pt,px,p>c ^ 1, in which case the gauge coupling constants can 
be disregarded, the fixed points px/pt = 1, Px/pt = and p^pt = 1, p\/pt = cease 
to be stable. Instead of them, the stable fixed point Rx = 3/4, = 3/8 [34j appears 
in the {Rx,R^) plane, where Rx = P\/pt and i?^ = p^/ Pt- In order to investigate the 
behaviour of the solutions to the renormalization group equations within the NMSSM, 
it is necessary to linearise the set of these equations in its vicinity and set = 0. As a 
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result, we obtain 



4 \^2^-^'^ v^"^ 8v^ 

+ I TT-KaO 7=-fL>£0 





Ai 



(7) 



where Rxo = -Ra(O), -R^ = R^{0), pto = Pt(0), Ai = — - — , A2 = — - — , and 

Ptit) = — ^ — • From (171), it follows that the fixed point Rx = 3/4, R^ = 3/8 arises 
1 + 7ptot 

as the result of intersection of two fixed lines in the {Rx,R^) plane. The solutions 

are attracted most strongly to the line -Rx H — ■i=R>c = „ ( 1 H — p V since Ai ^ A2. 

This line passes through three fixed points in the {Rx,R:^c) plane: (1,0), (3/4,3/8), and 
(0, 1). In the regime of strong Yukawa coupling, the fixed line that corresponds, in the 
{pt,p^,px) space, to the line mentioned immediately above is that which lies on the 
invariant surface containing a stable infrared fixed point. The line of intersection of the 
Hill and the invariant surface can be obtained by mapping this fixed line into the quasi- 
fixed surface with the aid of the set of renormalization group equations. For the boundary 
conditions, one must than use the values A^(0), x^(0), and h^{0) ^ 1 belonging to the 
aforementioned fixed line. 

In just the same way as infrared fixed lines, the infrared fixed surface is invariant under 
renormalization group transformations. In the evolution process, solutions to the set 
of renormalization group equations within the NMSSM are attracted to this surface. 
If boundary conditions are specified n the fixed surface, the ensuing evolution of the 
coupling constants proceeds within this surface. To add further details, we not that, 
near the surface being studied and on it, the solutions are attracted to the invariant line 
connecting the stable fixed point {px/ Pt = 3/4, p^^/ pt = 3/8) in the regime of strong 
Yukawa coupling with the stable infrared fixed point within the NMSSM. In the limit 
Pt, P>i, Pa ^ 1, the equation for this line has the form 
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(8) 



As one approaches the infrared fixed point, the quantities px and tend to zero: 
Pa ~ (pt — 7/18)^^/"*^^ and p^^ ~ {pt — 7/18)^/'^. This line intersects the quasi-fixed 
surface at the point 

P(to) = 0.82, P(to) = 0.087, pf P(to) = 0.178. 

Since all solutions are concentrated in the vicinity of this point for 
Yt{0),Yx{0),Y^{0) — ^ 00, it should be considered as a quasi-fixed point for the set 
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of renormalization group equations within the NMSSM. We note, however, that the 
solutions are attracted to the invariant hne (jH)) and to the quasi-fixed fine on the Hill 
surface. This conclusion can be drawn from the an analysis of the behaviour of the 
solutions near the fixed point {Rx = 3/4, = 3/8) (see ((7j)). Once the solutions 

have approached the invariant line -R\ H = - f 1 H ^ I, their evolution is 

governed by the expression (e(t))°-°^^, where e{t) = ptit)/pto- This means that the 
solutions begin to be attracted to the quasi-fixed point and to the invariant line (jH)) 
with a sizable strength only when Yi{0) reaches a value of 10^, at which perturbation 
theory is obviously inapplicable. Thus, it is not the infrared quasi-fixed point but the 
quasi-fixed line on the Hill surface (see Fig. 4) that, within the NMSSM, plays a key 
role in analysing the behaviour of the solutions to the renormalization group equations 
in the regime of strong Yukawa coupling, where all Yi{0) are much greater than cxq. 



3 Renormalization of the soft SUSY breaking pa- 
rameters 

If the evolution of gauge and Yukawa coupling constants is known, the remaining subset 
of renormalization group equations within the MNSSM can be treated as a set of linear 
differential equations for the parameters of a soft breakdown of supersymmetry. For 
universal boundary conditions, a general solution for the trilinear coupling constants 
Ai{t) and for the masses of scalar fields m'^it) has the form 

A{t)=eiit)A + mM,/2, (9) 
m^(t) = aiit)ml + bi{t)M^/2 + Ci{t)AMi/2 + di{t)A^ . (10) 

The functions ej(t), fi{t), ai(t), bi{t), Ci{t), and di{t), which determine the evolution of 
Ai{t) and mf{t), remain unknown, since an analytic solution to the full set of renormaliza- 
tion group equations within the NMSSM is unavailable. These functions greatly depend 
on the choice of values for the Yukawa coupling constants at the Grand Unification scale 
Mx- At the electroweak scale t = to, relations (jHI) and (fTn|) specify the parameters Af(to) 
and mf{to) of a soft breaking of supersymmetry as functions of their initial values at the 
Grand Unification scale. 

The results of our numerical analysis indicate that, with increasing Yi{0), where 

^ (ivrp' ^^^^"^ ^ (ivrp' ^"^^ ^""^^"^ ^ (ivrp' functions ei{to), ^(to), and di{to) 
decrease and tend to zero in the limit Yi{0) oo, relations and ()10j) becoming much 
simpler in this limit. Instead of the squares of the scalar particle masses, it is convenient 
to consider their linear combinations 

ml{t)=ml{t)+ml{t)+ml{t), (11) 
mlit) = 3mJ(t) 

in analysing the set of renormalization group equations. In the case of universal boundary 
conditions, the solutions to the differential equations for 9?t^(t) can be represented in the 
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same form as the solutions for mf (t) (see (fTO|l ): that is 



m^iit) = 3d^{t)ml + hi{t)Ml/^ + Ci{t)AMi/2 + di{t)A^. (12) 

Since the homogeneous equations for Ai{t) and 2ti(t) have the same form, the functions 
ai{t) and ej(t) coincide; in the hmit of strong Yukawa couphng, the rn^ dependence dis- 
appears in the combinations (fTT|) of the scalar particle masses as the solutions to the 
renormalization group equations for the Yukawa coupling constants approach quasi-fixed 
points. This behaviour of the solutions implies that Ai{t) and 9?t^(t) corresponding to 
^i(O) ^ «i(0) also approach quasi-fixed points. As we see in the previous section, two 
quasi-fixed points of the renormalization group equations within the NMSSM are of great- 
est interest from the physical point of view. Of these, one corresponds to the boundary 
conditions 1^(0) = 1a (0) ^ ai(0) and F^(0) = for the Yukawa coupling constants. The 
fixed points calculated for the parameters of a soft breaking of supersymmetry by using 
these values of the Yukawa coupling constants are 



p5fP(to)^1.77, pg,^P(to) ^ 6.09, 

p2r(^o)^-0.42, p«^P(to) ^ -2.^ 



(13) 



'A 



where PAiif) = A(^)/^i/2 and Pm^it) = Wlf/M^^^. Since the coupling constant x for 
the self-interaction of neutral scalar fields is small in the case being considered, Aj^{t) 
and 93T^(t) do not approach the quasi-fixed point. Nonetheless, the spectrum of SUSY 
particles is virtually independent of the trilinear coupling constant A^ since x — 0. 

In just the same way, one can determine the position of the other quasi-fixed point for 
Ai{t) and Tl]{t), that which corresponds to -Rao = 3/4, = 3/8. The results are 

pQfP(to) ^ 1.73, p^f (to) ^ -0.43, p^f (to) ^ 0.033, 

pQ^^(to) ^ 6.02, pg|P(to) ^ -2.34, pZ^'ito) ^ 0.29, ^^^^ 



where i?Ao = X\(0)/yi(0) and = Y^{0)/Yt{0). It should be noted that, in the vicinities 
of quasi-fixed points, we have p^^(to) < 0. Negative values of 93T|(to) lead to a negative 

value of the parameter m^to) in the potential of interaction of Higgs fields. In other 
words, an elegant mechanism that is responsible for a radiative violation of SU{2) ® U (1) 
symmetry and which does not require introducing tachyons in the spectrum of the theory 
from the outset survives in the regime of strong Yukawa coupling within the NMSSM. 
This mechanism of gauge symmetry breaking was first discussed in [SSj by considering 
the example of the minimal SUSY model. 

By using the fact that 2?t^(t) as determined for the case of universal boundary conditions 
is virtually independent of mg, we can predict aj(to) values near the quasi-fixed points 
(see 031 )• The results are 

1) -Rao = I5 R><o = 0, 

1 4 5 

ay{to) = au{to) = -, ai(to) = ag{to) = -, 02(^0) = ; 

2) -Rao = 3/4, R^ = 3/8, 

2 5 1 

ay{to) = 0, ai(to) = -02(^0) = g, s(^o) = g, an(to) = g • 



(15) 
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To do this, it was necessary to consider specific combinations of tlie scalar particle masses, 
such as m^— 2mQ, mQ+mfj—m'^+rnl, andm^— 2mf (at x = 0), that are not renormalized 
by Yukawa interactions. As a result, the dependence of the above combinations of the 
scalar particle masses on ttiq at the electroweak scale is identical to that at the Grand 
Unification scale. The predictions in (jlSj) agree fairly well with the results of numerical 
calculations. 

Let us now consider the case of nonuniversal boundary conditions for the soft SUSY 
breaking parameters. The results of our numerical analysis, which are illustrated in Figs. 
6 and 7, indicate that, in the vicinity of the infrared fixed point at = 0, solutions to 
the renormalization group equations at the electroweak scale are concentrated near some 
straight lines for the case where the simulation was performed by using boundary condi- 
tions uniformly distributed in the {At, Ax) and the (Tlf, plane. The strength with 
which these solutions are attracted to them grows with increasing Yi[0). The equations 
for the lines being considered can be obtained by fitting the numerical results displayed 
in Figs. 6 and 7. This yields 



For Y^{0) ^ ao solutions to the renormalization group equations are grouped near planes 
in the space of the parameters of a soft breaking of supersymmetry {At, Ax, A^) and 



It can be seen from Figs. 8 and 9 that, as the values of the Yukawa coupling constants 
at the Grand Unification scale are increased, the areas of the surfaces near which the 
solutions Ai{t) and OKf(t) are concentrated shrink in one of the directions, with the result 
that, at Yi{0) ~ 1, the solutions to the renormalization group equations are attracted to 
one of the straight lines belonging to these surfaces. 

The numerical calculations also showed that, with increasing Yi{0), only in the regime of 
infrared quasi-fixed points (that is, at -Rao = 1; R>cO = or at Rxo = 3/4, = 3/8) 
ei{to) and di{to) decrease quite fast, in proportion to 1/1^(0). Otherwise, the dependence 
on A and ttiq disappears much more slowly with increasing values of the Yukawa coupling 
constants at the Grand Unification scale - specifically, in proportion to (^^(O))"'', where 
S < 1 (for example, 6 = 0.35 — 0.40 at x = 0). In the case of nonuniversal boundary 
conditions, only when solutions to the renormalization group equations approach quasi- 
fixed points are these solutions attracted to the fixed lines and surfaces in the space of 
the parameters of a soft breaking of supersymmetry, and in the limit ^^(0) oo, the 
parameters Ai{t) and 971^ (t) cease to be dependent on the boundary conditions. 

For the solutions of the renormalization group equations for the soft SUSY breaking 
parameters near the electroweak scale in the strong Yukawa coupling regime one can 
construct an expansion in powers of the small parameter et{t) = Yt{t)/Yt{0): 



At + 0A47Ax = I.7OM1/2, 
mi + 0.147971^ = 5.76Mi% 



(16) 



{mlml,ml) (see Figs. 8-10): 



At + 0A28Ax + 0.022A^ = I.68M1/2, 
Tft + 0.128071^ + 0.022271^ = 5.77Mi^ 



(17) 




(18) 
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where and (3i are constants of integration that can be expressed in terms of v4j(0) and 
971^(0). The functions Vij{t) are weakly dependent on the Yukawa couphng constants at 
the scale Mx, and fij(O) = 1. They appear upon renormalizing the parameters of a soft 
breaking of supersymmetry from q ~ 10^^ — 10^^ GeV to g ~ rrit. In equations p8|) . we 
have omitted terms proportional to M1/2, ^1/2; ^i(0)Mi/2, and Aj(0)y4j(0). 

At X = 0, we have two eigenvalues and two corresponding eigenvectors: 

^=(3/7)' "=0 -3 

whose components specify {At, Ax) and (971^,971^). With increasing YtiQi) ^ ^a(O), the 
dependence on cto and (5q becomes weaker and the solutions at t = to are concentrated near 
the straight lines (^^(ai), A\[ai)) and (9Jlj 9)T|(/5i)). In order to obtain the equations 
for these straight lines, it is necessary to set A\{Q) = —3At{0) and 9Jt|(0) = —3971^(0) at 
the Grand Unification scale. At the electroweak scale, there then arise a relation between 
Aflto) and Ax{to) and a relation between 9Jt^(to) and 9)T|(to): 



A + 0.137Aa = 1.70Mi/2, ^^^^ 

(19) 

X — 0.1 yjivi^i^- 



971? + 0.137971^, = 5.76M2 



These relations agree well with the equations deduced for the straight lines at l^i(O) ~ 1 
by fitting the results of the numerical calculations ()16|). 

When the Yukawa coupling constant x is nonzero, we have three eigenvalues and three 
corresponding eigenvectors: 



1 



24 




whose components specify {At^Ax^A^ and (971?, 9)T|, 97t^). An increase in 
3 

yx(0) ~ 2Y'^(0) ~ ^^t(0) leads to the following: first, the dependence of Ai{t) and 

971? (t) on ao and /^o disappears, which leads to the emergence of planes in the space 
spanned by the parameters of a soft breaking of supersymmetry: 

A + 0.103Aa + 0.0124A^ = 1.69Mi/2, 
971? + 0.103971^ + 0.0124971?, = 5.78M?/2. ^ ' 

After that, the dependence on a\ and /3i becomes weaker at l^i(O) ~ 1. This means 
that, with increasing initial values of the Yukawa coupling constants, solutions to the 
renormalization group equations are grouped near some straight lines and we can indeed 
see precisely this pattern in Figs. 8-10. All equations presented here for the straight lines 
and planes in the 971? space were obtained at v4j(0) = 0. 

From relations (fTIJj) and (plj) . it follows that Atito) and 971?(to) are virtually independent 
of the initial conditions; that is, the straight lines and planes are orthogonal to the At 
and 971? axes. On the other hand, the v4^(to) and 971?,(to) values that correspond to the 
Yukawa self-interaction constant for the neutral fields are fully determined by the 
boundary conditions for the parameters of a soft breaking of supersymmetry. For this 
reason, the planes in the (Aj, ^4^, A^^ and (971?, 971^, 971?,) spaces are virtually parallel to 
the A„ and 9711 axes. 
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4 Conclusions 



In the strong Yukawa coupling regime in the NMSSM, solutions to the renormalisation 
group equations for Yiit) are attracted to quasi-fixed lines and surfaces in the space of 
Yukawa coupling constants and specific combinations of pi{t) are virtually independent 
of their initial values at the Grand Unification scale. For Yi{0) — > oo, all solutions to 
the renormalisation group equations are concentrated near quasi-fixed points. These 
points emerge as the result of intersection of Hill lines or surfaces with the invariant line 
that connects the stable fixed point for Yi ^ cij with the stable infrared fixed point. 
For the renormalisation group equations within the NMSSM, we have listed all the most 
important invariant lines and surfaces and studied their asymptotic behaviour for Yi ^ 
and in the vicinity of the infrared fixed point. 

With increasing Yi{0), the solutions in question approach quasi-fixed points quite slowly; 
that is, the deviation is proportional to {et(t)Y, where et(t) = Yt{t)/Yt{0) and 6 is calcu- 
lated by analysing the set of the renormalisation group equations in the regime of strong 
Yukawa coupling. As a rule, 6 is positive and much less than unity. By way of example, 
we indicate that, in the case where all three Yukawa coupling constants differ from zero, 
5 ~ 0.085. Of greatest importance in analysing the behaviour of solutions to the renor- 
mahsation group equations within the NMSSM at Yt{0),Yx{0),Y^{0) ~ 1 is therefore 
not the infrared quasi-fixed point but the line lying on the Hill surface and emerging as 
the intersection of the Hill and invariant surface. This line can be obtained by mapping 
the fixed points (1,0), (3/4,3/8), and (0,1) in the {R\,R^) plane for Yi ^ into the 
quasi-fixed surface by means of renormalisation group equations. 

While Yi{t) approach quasi-fixed points, the corresponding solutions for the trilinear 
coupling constants Ai{t) characterising scalar fields and for the combinations (t) of the 
scalar particle masses (see (jllj) ) cease to be dependent on their initial values at the scale 
Mx and, in the limit l^i(O) — > cxd, also approach the fixed points in the space spanned by 
the parameters of a soft breaking of supersymmetry. Since the set of differential equations 
for Aiit) and mf{t) is linear, the A, M1/2, and dependence of the parameters of a 
soft breaking of supersymmetry at the electroweak scale can be explicitly obtained for 
universal boundary conditions. It turns out that, near the quasi-fixed points, all Ai{t) 
and all ^f(t) are proportional to M1/2 and M'^^^, respectively. Thus, we have shown 
that, in the parameter space region considered here, the solutions to the renormalization 
group equations for the trilinear coupling constants and for some combinations of the 
scalar particle masses are focused in a narrow interval within the infrared region. Since 
the neutral scalar field Y is not renormalized by gauge interactions, A^{t) and ^^^^(t) are 
concentrated near zero; therefore they are still dependent on the initial conditions. The 
parameters At{to) and 93T^(to) show the weakest dependence on A and ml because these 
parameters are renormalized by strong interactions. By considering that the quantities 
971^ (to) are virtually independent of the boundary conditions, we have calculated, near 
the quasi-fixed points, the values of the scalar particle masses at the electroweak scale. 

In the general case of nonuniversal boundary conditions, the solutions to the renormal- 
ization group equations within the NMSSM for Ai{t) and 3Jl^(t) are grouped near some 
straight lines and planes in the space spanned by the parameters of a soft breaking of 
supersymmetry. Moving along these lines and surfaces as Yi{0) is increased, the trilinear 
coupling constants and the above combinations of the scalar particle masses approach 
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quasi-fixed points. However, the dependence of these couplings on Aj(0) and 3Jtf (0) dies 
out quite slowly, in proportion to {et{t))^, where A is a small positive number; as a rule, 
A ^ 1. For example, A = 3/7 at = and A ^ 0.0085 at ^ 0. The above is invahd 
only for the solutions Ai{t) and 931^ (t) that correspond to universal boundary condi- 
tions for the parameters of a soft breaking of supersymmetry and to the initial values of 
Rxo = 1, -RjiO = and Rxo = 3/4, = 3/8 for the Yukawa coupling constants at the 
Grand Unification scale. They correspond to quasi-fixed points of the renormalization 
group equations for Yi{t). As the Yukawa coupling constants are increased, such solutions 
are attracted to infrared quasi-fixed points in proportion to et{t). 

Straight lines in the {At, Ax, A^) and (971^, 9?l^) spaces play a key role in the analysis 
of the behaviour of solutions for Ai(t) and (t) in the case where Yt{0), 1a(0), Y^(0) ~ 1. 
In the space spanned by the parameters of a soft breaking of supersymmetry, these 
straight lines lie in the planes near which Ai{t) and OJl^(t) are grouped in the regime of 
strong Yukawa coupling at the electroweak scale. The straight hues and planes that were 
obtained by fitting the results of numerical calculations arc nearly orthogonal to the At 
and m'f axes. This is because the constants At{to) and 37t^(to) are virtually independent 
of the initial conditions at the scale Mx- On the other hand, the parameters A^{tQ) and 
9Jl^(to) are determined, to a considerable extent, by the boundary conditions at the scale 
Mx- At Rxo = 3/4 and i?^ = 3/8, the planes in the {At,Ax,A^) and (9Jlf , SJt^, 97t^) 
spaces are therefore parallel to the A^ and 3Jl^ axes. 
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Figure captions 



Fig. 1. The values of the Yukawa couphngs at the electroweak scale corresponding to the 
initial values at the GUT scale uniformly distributed in a square 2 < hf{0), A^(0) < 10. 
The thick and thin curves represent, respectively, the invariant and the Hill line. The 
dashed line is a fit of the values {pt{to) : Px{to)) for 20 < /i|(0), X^iO) < 100. 

Fig. 2. The values of the Yukawa couplings at the electroweak scale corresponding to the 
initial values at the GUT scale uniformly distributed in a square 20 < h^{Q), A^(0) < 100. 
The dashed line is a fit of the values (pt(io), Pa(^o)) for 20 < /i|(0), X^iO) < 100. 

Fig. 3. Evolution of the combination pt{t) + 0.506px{t) of the Yukawa couphngs from 
the GUT scale (t = 0) to the electroweak scale {t — to) for = and for various initial 
values /i2(o) - Fig. 3a, A2(0) - Fig. 3b. 

Fig. 4. Quasi-fixed surface in the Pa) space. The shaded part of the surface 

represents the region near which the solutions corresponding to the initial values 
2 < /i2(o),x2(0),A2(0) < 10 - Fig. 4a, 20 < /^^(o), ^^(O), A2(0) < lOO - Fig. 4b are 
concentrated. 

Fig. 5. Evolution of the combination pt + 0.720pA + 0.3330p;Kr of the Yukawa couphngs 
from the GUT scale (t = 0) to the electroweak scale (t — to) for various initial values 
/i2(0) - Fig. 5a, A2(0) - Fig. 5b, x'^{0) - Fig. 5c. 

Fig. 6. The values of the trilinear couplings At and A\ at the electroweak scale 
corresponding to the initial values uniformly distributed in the {At, Ax) plane, calculated 
at = and h'^{0) = X^{0) = 20. The straight line is a fit of the values (A(to), ^a(^o))- 

Fig. 7. The values of the combinations of masses 971^ and 9Jt^ at the electroweak scale 
corresponding to the initial values uniformly distributed in the {DJl^/My^,'Dyi1/My2) 
plane, calculated at = 0, h'^{0) = X'^{0) = 20, and ^(0) = Ax{0) = 0. The straight 
line is a fit of the values (3Jl?(io), ^l{to))- 

Fig. 8. Planes in the parameter spaces {At/Mi/2, Ax/Mi/2, A^/Mi/2) - Fig. 8a, and 
{m^/M^^^,ml/M^^^,ml/M^^^) - Fig. Sb. The shaded parts of the planes correspond 
to the regions near which the solutions at h^{0) = 16, A^(0) = 12, and x^(0) = 6 arc 
concentrated. The initial values Aj(0) and 971^(0) vary in the ranges —M1/2 < A < M1/2 
and < 97tf (0) < SM^^^, respectively. 

Fig. 9. Planes in the parameter spaces {At/Mi/2:Ax/Mi/2,A^/Mi/2) - Fig. 9a, and 
(3JttV-^i/2>^A/-^i/2'^y-^i/2) - Fig. 9b. The shaded parts of the planes correspond 
to the regions near which the solutions at hf{0) ~ 32, A^(0) = 24, and x^(0) = 12 are 
concentrated. The initial values Aj(0) and M|(0) vary in the ranges —M1/2 < A < M1/2 
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and < 071^(0) < ^Ml,^, respectively. 



Fig. 10. Set of points in planes (0.0223(A^/Mi/2) + 0.1278(Aa/Mi/2), A/M1/2) - Fig. 
10a, and {^m2?,{ml/Mlf^) + 0.1278(9Jr^/M2/2), T^/M^^^) - Fig. 10b, correspond- 
ing to the values of parameters of soft SUSY breaking for hf{0) = 32, A^(0) = 24, 
x^(0) = 12, and for a uniform distribution of the boundary conditions in the parameter 
spaces {At,Ax,A^) and (271^^ ajt|, 071^). The initial values ^^(0) and 971^(0) vary in the 
ranges — M1/2 < A < M1/2 and < 971^(0) < 3My2, respectively. The straight lines in 
Figs. 10a and 10b correspond to the planes in Figs. 9a and 9b, respectively. 
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Fig.9b. 
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Fig.lOa. 
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